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' Abstract 

, ^ ' Chains of first-order SUSY transformations for the spin equation are studied in detail. It 

. is shown that the transformation chains are related with a polynomial pseudo-supersymmetry 

\ of the system. Simple determinant formulas for the final Hamiltonian of a chain and for so- 

lutions of the spin equation are derived. Applications are intended for a two-level atom in an 
. electromagnetic field with a possible time-dependence of the field frequency. For a specific form 

ly-^ \ of this dependence, the time oscillations of the probability to populate the excited level disap- 

0^ ■ pear. Under certain conditions this probability becomes a function tending monotonously to a 

constant value which can exceed 1/2. 
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1 Introduction 



It is well-known that complex quantum systems with a discrete energy spectrum (i.e., an atom in 
\ an external electric field) can be placed in a special dynamical configuration in which only two 
qh' stationary states are essential. In those cases the description of a system's evolution does not require 
a consideration of the entire Hilbert space. A good approximation is achieved in case the whole 
Hilbert space is replaced by a two-dimensional linear space. Systems of this kind are called two-level 
systems, whose evolution is governed by a special system of two differential equations, known as the 
"spin equation" . 

The spin equation arises in many areas of theoretical physics and thus finds a wide range of 
applications, for instance, the semi-classical theory of laser beams jj, the absorption resonance and 
nuclear induction experiments j2], the behavior of a molecule in a cavity immersed in electric or 
magnetic fields rSj. Particularly, we would like to stress the importance of the spin equation in view 
of its possible applications to quantum computations 0. 

The first exact solution of the spin equation was obtained by Rabi |5j . This solution found an 
extensive physical application in the study of various properties of two-level systems. This fact shows 
a great importance of exactly solvable models involving the spin equation. In a recent paper jH] , some 
general properties of the spin equation were studied. The authors of jB] give an overview of the known 
results and obtain new classes of exact solutions to this equation. Among the methods of analysis 
of exactly solvable spin equations, the method of intertwining operators 13^0 plays a special role, 
since it reveals a new type of symmetry related with two-level systems 13111, called the polynomial 
pseudo-supersymmetry. 

In this respect, we would like to remind that the ideas and methods of supersymmetry (SUSY) 
are recently coming into use in relativistic and nonrelativistic quantum-mechanical problems. Su- 
persymmetric quantum mechanics (SUSY QM), originated from the simplest quantum field model 
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in connection with the problem of a spontaneous SUSY breaking jTUj, finds more and more extensive 
apphcations in theoretical physics (see, e.g., [H]). It was discovered that diverse aspects of SUSY 
QM are closely related with the method of Darboux transformations, well-known in the soliton the- 
ory Numerous works are devoted to the relation between SUSY QM and the technique of 
intertwining operators [ini EI • 

This article is a continuation of a previous works devoted to the study of exact solutions of two- 
level systems and the peculiarities related to the corresponding underlying supersymmetry [ZI-|ni. 
Our main objective is to show that a repeated use of simple transformations introduced in the above 
papers gives rise to a wide range of exactly solvable interactions for the spin equation. It turns out 
that such a study is convenient in the analysis of equations describing two-level systems written in 
the form of a one-dimensional stationary Dirac equation with a non-Hermitian Hamiltonian of a 
special form, in which time plays the role of a spatial variable [7]. In this respect, we can mention 
that non-Hermitian interaction arises, in particular, in various problems of field theory, statistical 
mechanics (see, e.g., fj^) and nuclear physics. In connection with the discovery of a large class 
of complex potentials possessing the so-called PT-symmetry (see, e.g., |16j), which often leads to 
Hamiltonians with a purely real spectrum, recently an attempt was made once again to construct 
a complex extension of quantum mechanics ^7]. In this respect, it is worth mentioning the paper 
|18j . where the authors show that the use of non-Hermitian (so-called quasi-Hermitian) operators 
does not contradict the basic principles of quantum mechanics. It should also be noted that such a 
generalization was not started from scratch, and, in fact, partially goes back to classical works [T9] . 

In this respect, the role of supersymmetry in quantum mechanics may increase |H], since, as 
was recently shown, SUSY transformations allow one to remove singular points from continuous 
spectra of non-Hermitian Hamiltonians and to convert non-diagonalizable Hamiltonians into 
diagonalizable ones [21]. This may enlarge the set of admissible super-Hamiltonians, thus giving an 
opportunity to extend the theory by new types of interaction. 

In it was shown how to construct matrix-differential intertwining operators preserving the 
special form of a non-Hermitian Dirac Hamiltonian, and in ,9] some results were announced related 
to transformation chains. The technique of simple (first-order in derivatives) intertwining operators 
adapted to the special case of Rabi oscillations made it possible to discover an interesting physical 
effect Namely, as distinct form a constant field frequency when the probability to populate the 
excited level oscillates with time (Rabi oscillations), there exist such types of time dependence of 
the field frequency that this probability ceases to oscillate and becomes a monotonously increasing 
function of time tending to a value which may exceed 1/2 (in particular, it may be equal to 3/4 [S]). 
This property gives the hope that inverse population may be observed in an ensemble of two-level 
atoms placed in such type of field, and they may exhibit lasing properties. This result was obtained 
when a single Darboux transformation was applied to a two-level atom. 

One of the aims of the present work is a detailed analysis of the above-mentioned transformation 
chains. In order to make this work self-contained, we show first of all (see the following section) how 
the Schrodinger equation for a two-level atom interacting in the rotating wave approximation with 
the electric component of an electromagnetic field can be reduced to a one- dimensional stationary 
Dirac equation with an effective non-Hermitian Hamiltonian where time plays the role of a spatial 
variable. This is accomplished for a general time-dependence of the field frequency. 

The third section is devoted to supersymmetric constructions. It is shown that the transformation 
chains introduced in [7| |H] are related with a polynomial pseudo-supersymmetry of the system. This 
fact indicates the presence of this type of symmetry in a two-level atom. 

The introductory part of the fourth section contains a brief review of previous results which are 
required in the following sections. It is shown that in the case of various factorization constants 
the use of formulas similar to the Crum-Krein formulas |22| for the Schrodinger equation permits 
us to express the potential, obtained as a result of an n-fold SUSY transformation, through n-order 
determinants of the transformation functions, thus avoiding the use of recurrent formulas, which 
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require the knowledge of solutions at all the intermediary steps. We obtain the form of transformation 
functions that preserves both the specific matrix structure of the initial potential and its real- valued 
character. We also study the case of coinciding factorization constants and obtain a sufficiently 
simple realization of such transformation chains. 

As a numerical application of our general scheme, in the final (fifth) section we examine a two- 
fold transformation for a two-level atom in an external electromagnetic field. We discover that in 
this case the disappearance of oscillations in the time-dependence of the probability to populate the 
excited level may also occur, i.e., at certain conditions this probability may acquire a monotonous 
time-dependence. It is established that, in contrast to the previously published results |HI, this 
effect emerges for two kinds of behavior of the detuning of the external field frequency from the 
resonance value. In these cases, the probability values to populate the excited level after a two-fold 
transformation coincide with either the maximal or the minimal values of the probability oscillations 
obtained ClS Sb result of a one-fold transformation. 



2 Two-level atom in external field 

Since the results obtained in this paper are mainly applied to a two-level atom in an external elec- 
tromagnetic field, in this section we briefly review some properties of this system which are required 
in the following sections. 

In the non-relativistic approximation, the evolution of a state vector of a two-level atom in an 
external electromagnetic field with the electric component E (t) = eE^ cos [uj [t) t], where e is a unit 
vector of the field polarization, is described by the Schrodinger equation 

th^{r,t) = {Ho + H,)^{r,t) . (1) 

Here, Hq is the conventional Hamiltonian of an atom in the absence of an external field, Hoipnlr) = 
£„^/'„(r), n = 1,2] Hi describes the interaction of the atom with the electric field component in the 
dipole approximation. 

Hi = \q\ Eq (e ■ r) cos [u {t) t] ; 

g < is the charge of electron; the dot over a symbol stands for differentiation with respect to time. 
We assume that the field frequency, ui = ui {t), and thus also the detuning of the frequency from the 
resonance value, to'21 = c<j2 — cui, 

5{t)=uj{t)-U2i, (2) 

are functions of time. Decomposing the solution ip (r, t) of equation ((T)) in the basis ipn (r) exp (— 

= ^n/T^i n = 1,2, we obtain the well-known [2^ set of equations for the decomposition co- 
efficients Cn, n = 1,2. Introducing the notation Ai^2 = exp (=F*f^) Ci,2 and neglecting the term 
exp [i{uj2i + uj(t))t], rapidly oscillating in comparison with exp [i6t] (the so-called rotating wave ap- 
proximation I24j ) ■ we arrive at the following set of equations: 

lAi - fAi = , iA2 + fA2 = ^Ai . (3) 

Here, ^ = ^Eod2i, 

f = f(t) = ~mt)], (4) 

and d2i is a matrix element of the dipole momentum operator, di2 = 621 = \ (e • r) \ip2)- If 

does not depend on time (consequently, f = = const), then solving equations © with the initial 
conditions A2 = and Ai = 1 for t = we determine the probability to detect the system at the 
moment t in the excited state on condition that at the initial time moment t = the system was at 
the ground state: 

\A2it)\' = ^[l-cosi2nt)], n' = f + e. (5) 
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This result is well-known in quantum optics j2Sl I21|- Formula (0) describes the so-called Rabi 
oscillations, and the value equal to 2^, is known as the Rabi frequency. 
Introducing the notation 

Voit) = ta2fo (t) , (6) 

7 = iai, E = ^, = {Ai,A2y (the symbol^ stands for transposition, while 0"i 2,3 are the standard 
Pauli matrices) we rewrite the system of equations © at / = /o in the matrix form 

/lo^ = , ho=7j^ + Vo (t) . (7) 

Equation (|7j) has the form of a one-dimensional stationary Dirac system where time plays the role 
of a spatial variable and Vq [t) is a matrix- valued potential determined by the function /g (t), which 
we shall call the potential. Following the well-established terminology, we call ho the Hamiltonian 
even though in the current case it does not correspond to any quantum-mechanical system. By 
construction, the parameters /o and E are real. 



3 Polynomial pseudo-supersymmetry of a two-level system 

Let us assume that the solutions of equation ((7j) are known and it is necessary to find solutions of 
this equation with another potential, 

/ii$ = E<l>, h,=^j^ + Viit) , (8) 

where Vi{t) = ia2fi (t). This problem can be solved by finding such an operator Lo,i that obeys the 
following operator equality (intertwining relation): 

^0,1^0 = hiLo^i . (9) 

In this case, solutions of equation (jHJ can be found by applying the operator Lq i to solutions of the 
initial equation ((Tj), $ = Lq i"^. 

The intertwining relation © provides the basis for the general concept of transformation operators 
(see, e.g., 123), and, in particular, operators of Darboux transformations [2n]- At the same time, 
the existence of the intertwining operator Lo,i allows one to construct a (generally polynomial) 
supersymmetry algebra [2^1 related to equation ((7j), thus revealing its internal supersymmetric nature. 

Usually in optical problems, as distinct from quantum-mechanical problems, the introduction of 
a Hilbert space is unnecessary. In this case, symmetry operators are defined on the space of solutions 
of the corresponding equation without introducing any inner product. In the case of a two- level 
atom, the equation in question is the differential equation (jTj). For supersymmetric constructions, 
we shall need the notion of operator conjugation, which will be introduced in a formal way. The 
operation of formal (Laplace) conjugation obeys the standard rules {ABy = B^A\ {d/dtY = —d/dt 
and corresponds to the transposition of a matrix accompanied by the complex conjugation of its 
elements. The operator A is called Hermitian if A"^ = A. It is easy to seen from ((Tj) that ho is 
non-Hermitian. Therefore, along with equation ((Tj) there also exists its adjoint form 

hl^ = E^. (10) 

Another relation that we shall need is obtained by the conjugation of the intertwining relation Q 

Ll,hl = hlLli- (11) 
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The operator Lq i, obviously transforms solutions of equation ((Tj) into solutions of equation (jSI), while 
the operator L^^ realizes the backward transformation, i.e., a transformation from solutions of the 
equation 

h\(^ = E$ , $ = 4 1^ 

into solutions of equation pOj) . We shall also assume that the functions \1/ and $, $, as well as 
the operators ho, and hi, h\ are interrelated by an operator J 

^ = , $ = J$ , 

which is assumed to have the properties 

= Jho^iJ , J2 = ±1 , Jt = ± J . (12) 

From these equations it follows that a superposition of conjugation and the operation J is a pseudo- 
conjugation, while the operators ho and hi are pseudo-Hermitian ^Tj . 
Relations (fTT|) and (fT^ yield the equahty 

JIjQiJhi = hoJZjQ^J . 

Thus, the operator JLq ^ J transforms solutions of the equation with the Hamiltonian hi into solutions 
corresponding to the Hamiltonian ho, while the superposition JLq^JLq^i transforms solutions of 
equation (|7j) into solutions of the same equation, and, therefore, it is a symmetry operator of this 
equation. In exactly the same manner, the operator Lq^iJLq^J is a symmetry operator of the 
equation with the Hamiltonian hi. Since both these operators are differential, one can expect that 
the indicated superpositions are polynomials of the corresponding Hamiltonians, a property which 
we shall demonstrate. 

It is easy to see that in the case of Hamiltonian (|7j) we have J = ax- In the authors constructed 
matrix-differential intertwining operators preserving the specific form (j?)) of the non-Hermitian Dirac 
Hamiltonian ho and expressed the symmetry operators JLq ^ J Lqi and LqiJLq^J of equations ^ 
and ((HI) in terms of the corresponding Hamiltonians: 

JL+i JLo,i = hl-Al, Lo,i JL+i J = /i? - . (13) 

The constant matrix Ai = diag(Ai, — Ai) in (fT!?|l is called the (matrix) factorization constant (Ai 
is also called the factorization constant). Formulas (fT^ present a generalization of the factoriza- 
tion properties of transformation operators that take place in the case of Hermitian one-component 
Hamiltonians ^28j. 

The presence of a polynomial pseudo-supersymmetry in the system under consideration is related 
with the possibility of a repeated use of the above one- fold transformation, i.e., with the realization of 
transformation chains. Let a sequence of transformation operators Lo,i, -Z^i,2, • • • , -^n-i,n intertwines 
non-Hermitian Hamiltonians hQ, hi, hn and all the operators h^ have potentials of the form 
Vfc if) = io'2fk (t)- In this case, every h^ has the property 

h+ = JhkJ, k = 0,...,n. (14) 

If we are interested only in the resulting action of the chain, we can introduce the operator Lo,„ = 
Ln-i^n ■ ■ ■ Li^2Lo,i- Due to the intertwining relations which hold for each separate operator L^-i^k 

Lk-i,khk-i = hkLk-i^k , k = 1, . . . ,n, 
the operator Lo,n obeys the relation 

(15) 
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Conjugating (fT^ and using (HH), we obtain 

JLlJK = h,JLlJ. (16) 
Moreover, for every operator there holds the following factorization: 

'JLk~i,k'J^k-i,k = hl^i - -^l, k = l n 



with Afc = diag(Afc, -Afc). 

Let us consider the superposition 

J^t,nJ^O,n = JLq iLi2 ■ ■ ■ L^-l,nJ ^n-l,n ■ ■ ■ -^1,2-^0,1 • 

A subsequent use of separate factorization relations of type (\Vi\ji and intertwining relations of type 
Q yields 

JL+„ JLo,„ = {hi - Al){hl -Al)... {hi - Al) . (17) 
In a similar way, we obtain the factorization of the same polynomial of the operator hn 

L,,nJLlJ = {hi - Al){hl -Al)... {hi - Al) . (18) 

Let us introduce the following 4x4 matrices: 

/lo \ ^ _ ^ 0\ ^ _ /O JL+ J 



^"lo h^) ' ^^"U," oj ' vo 

It is easy to see that the intertwining relations (fTH|) and (fTHIl are equivalent to the commutation 
relations 

\Q^,H] = [Q2,H] = Q, (19) 
while formulas p7p and psp can be rewritten in the form 

Q1Q2 + Q2Q1 = {H' - Tj) {H^ -Tl)... {H' - Tl) (20) 

where = diag(Afc,Afc). It is also obvious that the operators Qi and Q2 are nilpotent. Formulas 
(fT^ and (pUj) demonstrate that the operators H, Qi and Q2 close a polynomial superalgebra thus 
revealing a hidden polynomial supersymmetry of the system described by the Hamiltonian ho. For 
the first time, superalgebras related with non-Hermitian Hamiltonians were introduced in |2ZI. From 
this standpoint, a polynomial pseudo-supersymmetry can be associated with a two-level system. 



4 Chains of transformations 
4.1 General remarks 

Supersymmetric constructions of the previous section are valid if there exist transformation operators 
Lo,i,. . . ,Ln-i,n preserving the special form ((Zj) of a non-Hermitian Dirac Hamiltonian at each step of 
transformation where the potential matrix has the form In jH] the authors constructed one- fold 
matrix-differential intertwining operator Lq i having the necessary properties. In this section, we 
shall sketch the results of [8^ that we shall need in the following sections. 
Let us assume that we know n matrix-valued solutions of the equation 

{jdt + Vo)Uk=UkAk, k = l,...,n. (21) 
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This equation always has a solution, 



Uk Uk 
Uk ~Uk 



(22) 



corresponding to the factorization matrix = diag {iRk-, —iRk)- The functions Uk-, Uk have the same 
absolute value, and the sum of their phases is a constant, which we denote by C; see eqs. IH^ . (j211))- 
We choose the matrix-valued function Ui as a transformation function for the first transformation 
step. The functions ui, ui defining the matrix Ui are expressed in terms of a real function gi. 



Ul 



Ul 



iqi 



exp 



exp 



Ri 



Ri 



1-9? 



1 



dt 



dt 



(23) 



(24) 



which, in turn, should be found with the help of a real solution xi of the second-order equation 



Xi + 



Xi 



as follows: 



Ri 
fo 



fo 



Xi 



(25) 



(26) 



2/o' foXi ■ 

The transformation operator of the first transformation step Lq i and a new potential Vi have the 
form 



Vi = 



dt 



Di = Wi. 



(27) 
(28) 



The potential matrix Vi 



has the form ©, where /o should be replaced by /i. The function /i 



and the matrix Wi defining the operator Lq i are expressed in terms of the same function gi, 

4i?igi 



/i 



W 



1,2 ^0 ' 

l + 9i 

diag (1^1,^2) 



Wi 



-tfo + Ri 



W2 



(29) 
(30) 



If we choose Ai to be purely imaginary, the parameter Ri becomes real and the differential 
equation has real coefficients. Due to the fact that an equation with real coefficients always 
has real solutions, we conclude that for any real- valued function foit) one can obtain a real- valued 
function fi{t) thus realizing a transformation with a real- valued resulting potential. 

We now observe that the functions Vk = LoiUk, k = 2, . . . ,n are matrix- valued solutions of equa- 
tion (|2H) with potential (j28|) . At the next step of transformation, we choose V2 as a transformation 
function, thus obtaining the operator Li 2 and potential V2 



L - ^ 



Vo 



W2, W2 = K2 K2 

Vi + [^,W2] = Vo+b,D2] 



(31) 



D2 = Di + W2. 



The operator Lo,2 = -^1,2-^0,1 is a differential operator of second order with matrix-valued coefficients, 
intertwining the Hamiltonians and /i2 = idt + V2, -^^0,2^0 = ^2-^0,2- 
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Continuing this process, one can construct the n-th order operator 

(32) 

which transforms the Hamiltonian Hq into the Hamiltonian hn with the potential 

K = K-i + [7,^n] =H + [7,^n], Dn = Dn-i + Wn, Do = , (33) 

where 

Wn = yr.y-\ yn = ^0,n-lWn • (34) 

The derivation of the transformed potential Vn by the recurrent formulas (jHHj) and (jH^ involves 
a calculation of all the intermediary potentials and transformation functions, which requires a large 
amount of computational work. Below, we will show how to overcome this difficulty in the case when 
all the factorization constants . . . , Rn are different from each other. For coinciding factorization 
constants, we will show that equation may be solved in quadratures for any initial potential /o 
which also yields an essential simplification of the method. 



4.2 Transformations with coinciding factorization constants 

Let us consider a chain of two transformations. In order to obtain the potential /2, we have to make 
the replacement Ri ^ R2, Qi ^ (I2, fi ^ f2, fo ^ fi in (EHI)) and then fi in the obtained equation 
should be substituted by using the same equation ^2^ . This yields 

4Rigi 4R2g2 

In (j35|) the function q2 is expressed through the solution X2 of equation (j25|) with /o = /i and 
Ri = R2 by formula (jSH)), where the necessary replacements should be made. A generalization of 
this procedure to the case of a chain of n transformations is obvious. However, the form of both 
the potential and solutions to equation ()25|) after each transformation step will become more and 
more complicated. It is remarkable that, in case the transformation chain is realized by coinciding 
factorization constants, equation can be integrated explicitly. This makes it possible to get rid 
of this equation and express the new potential /2 in terms of the potential /i alone. 

Note, first of all, that according to for Ri = R2 and q2 = qi or q2 = qi^ , the potential f2 

assumes the value of the initial potential /o, i.e., the second transformation step is equivalent to a 
reverse transformation. The case q2 = qi corresponds to a trivial transformation, for which all of the 
three potentials remain the same, f2 = fi = fo, whereas for q2 = qi^ one obtains a nontrivial result. 
Therefore, for R2 = Ri the inverse transformation corresponds to q2 = qi ■ Then using formula 
dSni) for q2, in which we have to put qi = q2, fo = fi, Xi = X2, we obtain the following differential 
equation with respect to the function X2'- 



which can be easily integrated: 



-1 ^1 ^2 /oc\ 



X2 = Cexp {R,-^-^]dt. (37) 




Here, C is an integration constant. Equation with fixed values of /o and Ri, being a second-order 
differential equation, has two linearly independent solutions. Using the fact that the Wronskian of 
any two linearly independent solutions to equation is constant |l2n|, we can obtain a solution 
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of this equation, which is hnearly independent from (j37p . This permits us to find the potential /2 
different from /o, i.e., to accomphsh the second nontrivial step of transformation, with the same 
value of the factorization constant. To this end, let us define the function X2 as the second solution 
of equation satisfying the condition W (x2, X2) = 1, where W denotes the Wronskian. Then the 
function X2 is expressed in terms of X2 as follows: 



X2 = X2 



jx2^dt. (38) 



Thus, in order to obtain a real potential by a two-fold transformation with coinciding factorization 
constants, it is now unnecessary to solve equation at the second step of transformation, since 
the structure of its solution leading to a nontrivial result is determined and expressed in terms of 
quantities known as a result of the first transformation step. 

Let us express the result of consecutive application of transformations with coinciding factoriza- 
tion constants in the general form. A solution of equation ()25|) subject to the substitution /o fn-i 
and Ri — > Rn corresponding to the inverse transformation such that /„ = fn-2 has the form 

Xn = Cexp [ i Rn-i - /""^ - I dt , (39) 
J \ 2/„_i j 

where Rn = Rn-i- By analogy with we obtain the second solution Xn of this equation linearly 
independent from (jSHj), corresponding to another nontrivial step of a many-fold transformation. The 
new potential can be found by formula subject to the substitution /o fn-i and qi qn, 
where 

-Rn fn—1 Xn r An\ 

qn = -r iTTT- - • (40) 

Jn-l Jn-lXn 

Thus, the use of transformation chains with coinciding factorization constants allows one to easily 
obtain new exactly solvable potentials. 



4.3 Transformations with different factorization constants 



As we have already mentioned, the explicit solution of equation (^3)) becomes problematic in the case 
of a many-fold transformation. The integration of this equation can be avoided by using formulas 
similar to the Crum-Krein formulas [22j, well-known for the Schrodinger equation, which allow one 
to express the potential after an n-fold transformation in terms of determinants of transformation 
functions. This question was considered in detail for a matrix Schrodinger equation [201 ^^"^ ^ Dirac 
equation J3] with Hermitian Hamiltonians. Below, we will investigate the peculiarities arising during 
the use of these formulas for the Dirac equation ((7j) with the effective non-Hermitian Hamiltonian 
Hq and transformation functions of a special type, which do not change the form of a Hamiltonian 
after a transformation and obey conditions necessary for the transformed potential to remain a real 
value. 

We start by reminding some results concerning chains of transformations for the Dirac equation 



obtained in [1311201 • According to these articles, the entries d^^ 
are expressed in terms of the transformation functions Uk, k 



i,j = 1, 2 of the matrix D„ in ()33p 
n as follows: 



jn 



I Pij (Ml , ■ ■ ■ Mn 



\P{U 



1; 



Here, 
1,... 



P(Wi, 
n, 



,Un) is the Wronsky matrix constructed from the matrix- valued functions Uk, k 



( 



1) 



U2 
U2' 



(n-1) 



1) 



(41) 
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the matrix Pij {Ui, . . . , W„) has a structure similar to that of P (Ui, . . . , W„) with the exception that 
the final row consisting of matrices W^" ^'^ is replaced by matrices , = 1, . . . , n: 



Pij (Wl, • • • , Wn) 



K2 



(n-2) y{n-2) 



The matrices i,j = 1, 2 follow from by the replacement of the j-th row by the z-th row of 

the matrix U^"'^ . 

In jnni it was shown that the action of the operator Lo,n p2|l on the function = {iPie,4'2e) 
yields a vector = i<PiE, <p2E)'^ with the entries 



IJjE 



\PjE ■ ■ ■ Mn) 



In (1121) the matrix PjE (Ui, . . . ,Un) has the form 



PjE ■ ■ ■ Mn 



( Ui 



u 



U2 

(n-l) 



1,2 



lAn 



(42) 



^E 



u 



(n-l) ,T,(n-l) 



(43) 



where Ul is the j-th row of the matrix Uk- These formulas are obtained as a closure of the recur- 
sion scheme defined by equations ()33|) and ()34j). This algorithm will preserve the particular form 
of potential (0) provided that the transformation operators L^-i^k preserve the form (j^^ of the 
transformation function at every transformation step, a property which we will now demonstrate by 
induction. 

Let the transformation function Ui of the first transformation step have the required form 
We also assume that this property takes place at the {k — l)-st transformation step: 



V, = L 



Q,k~2Uj 



^3 ^3 



j = k-l,k,...,n. 



■^3 ^3 

Then the transformation function at the A;-th step is expressed as follows: 



yk 



^k-2,k-l 



^k-2,k-l 



(44) 



(45) 



Using (jUJ, we find, just as we expected, that imphes 

yk = 



Yk Yk 
Jk -Jk 



Vk-l 

Jk = Vk Vk . 

Vk~l 



It has been noted that for a purely imaginary parameters the condition necessary for the 
transformed potential to be a real value is the condition that the absolute values of the elements of 
the matrix-valued transformation function U be equal. One can see that the set of functions (j^^ 
actually meets this condition. Therefore, these functions allow one to realize a transformation chain 
which preserves both the special form of the potential and the condition for the function fk to be 
real- valued. 
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We are now going to show that for the choice of transformation functions in the form (|22|) the 
matrices of (2n)-th order P (Ui, . . . Mn), Pij iMi, ■ ■ ■ Mn) assume a block form with one of the blocks 
being a zero matrix. As a result, the determinants \P (Ui, . . . ,lAn) |, \Pij (pfi, ■ ■ ■ Mn) \ are reduced to 
the product of two determinants of ra-th order. Let us consider in more detail how this happens for 
n = 2, and then present the result for an arbitrary n. 

In what follows, the elements, rows and columns of a determinant are understood as the elements, 
rows and columns of the corresponding matrix. 

According to ()4H). the determinant \P {Ui, . . . ,Un) I for n = 2 has the form: 



\P{UiM2) 



Ui Ui U2 U2 
Ui -Ui U2 -U2 



U 



U 



1 



U2 

U'n 



U 



2 
U2' 



(46) 



Subtracting from the second column the first one, and from the forth column the third one, and 
transmuting the rows and columns, we reduce the determinant fl46|) to the product of two second- 
order determinants: 



\P{UiM2) 



It is clear that in general we have \P (Ui, . . . ,Un) 



Ui 


U2 












u[ 


u'2 










Ml U2 


Ui 


U2 


-2mi 


-2U2 




u'l u'2 


u[ 


u'2 


-2ui' 


-2S2' 







Ui U2 

u'l u'2 



-2)"' P1P2, where 



where 





Ul 


U2 


Un 




Ul U2 


Un 




< 


u'2 


■■ < 




u'l u'2 


... < 


Pi = 


uf'^ 


(n-2) 
U2 




, V2 = 


~(n-2) ~(n-2) 
Ul U2 


~(n-2) 
■ ■ ■ Un 




(n-1) 
Ml ' 


(n-l) 
u'2 


.. ut-"^ 




U\ U2 


~{n-l) 
. . . Mn 


at ions 
): 


allow 


one to simplify the determinants of the matrices Pu 




.Mn)\ = 


(-2)"rip2, |P22(Wi,...,W„)| = 


[-2Tpir2, 




Ul 


U2 


Un 




Ul U2 


Un 




u'l 


u'2 


■ W'n 




u'l u'2 


... M^ 


ri = 


uf''> 


{n-2) 
U2 


{n-2) 
. M„ 


, r2 = 


~(n-2) ~(n~2) 
Ul U2 


~(n-2) 
. . . Un 






(n) 
U2 


(n) 

. u\ 




~(n) ~(n) 
u\ ' U2 


~{n) 

... u\' 



The determinants of the matrices P12 (Ui, . . . ,Un) and P21 (^i, • • • Mn) for the transformation func- 
tions are equal to zero. Threfore, the matrix Dn from ()33|) becomes Dn = diag (^^) and Vn 
can be written in the form 

Vn = Vo + AVn, AVn = ta2Afn, A/„ = /„ - /o = - 



ri 
Pi 



P2, 



Thus, we have shown that the use of the determinant technique for an n-fold transformation of 
the spin equation, which is similar to the Crum-Krein approach to the Schrodinger equation, permits 
us to reduce the method to calculating only n-th order determinants. 
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5 Transformations of the Rabi oscillations 



The Rabi oscillations occur at /o = const. In !B| a detailed analysis of solutions to the set of equations 
((Tj) constructed with the help of a one-fold transformation is given. Here, we shall study the effect 
of a two-fold transformation. 

Let us first consider the probability to populate the excited level obtained as a result of a two-fold 
transformation with coinciding factorization constants R2 = Ri- For /o = const, equation takes 
a very simple form: 

Xi + zu\, = 0, w^ = fl-Rl. (47) 

Solutions of this equation have various properties depending on whether the value of cc^ is posi- 
tive, negative or zero. We have found |H] that the Rabi oscillations after a one-fold transformation 
disappear for tu = 0. In this case, the general solution of equation (P7j) is a linear function of 
time, Xi = + -B, which in the case A = 2Bfo yields the following expression for the function qi: 
gi = 1 — 2/ (2/ot + 1). The potential after the one-fold transformation reads 

/i(t) = /o 



1 + 4/o^t2 

Let us make the next transformation with R2 = Ri- The function q2 is found from at n = 2, 
where X2 is determined by formulas (|!?7j) - PH|) . Using we get the following expression for the 
transformed potential: 

r ^ /o(-648(-l + 4m - 243(1 + + {C - 6/ot(9 + 4^)') 

81(l + 4/2t2)2 + (C-6/ot(-3 + 4/2t2))2 ' ^ ) 

where C is a constant. 

Let us determine the form of the external electric field E = E{t) = Eq cos [tuj2 {t)] corresponding 
to the potential /2 {t) ()48|) . We may observe an interesting feature of the behavior of E{t) with 
time. From ^ and Q, it follows that the frequency 002 {t) is expressed through the function /2 (t) 
as follows: ^ 

U2 (t) = 0021 + S2 (t) , 52 (t) = 7 / /2 (r) dr . (49) 

t Jo 

In general, the field E{t) is an oscillating function with the amplitude Eq (see Fig. 1) except for 
a special time dependence of 62(t). In particular, it may happen that S2(t), as found from (jlH)), is 
expressed in terms of an inverse trigonometric function shifted by a constant cq. If cq + UJ21 = 0, the 
function cos[ta;2(t)] may become polynomial. For instance, substituting in (jlHj) C = and /o = —1, 
we obtain for uj2{t) 

4 fSt^-Qt 
uj2\t) = LJ21 — 2 H — arctan 



t \12t^ + 3 

Thus, for 1^21 = 2 the field E becomes the polynomial 

^ , , ^ f 5184(1 + 16t2 + 48t^) 288(2 + Tt^ - 4t^ 



(9 + 108^2 + 48t4 + 64t6)2 9 + 108^2 + 48^^ + 
plotted in Fig. 2. 

Another interesting effect has been observed for a one- fold transformation 8]. It consists in the 
disappearance of oscillations in the time-dependence of the excited level population. For a two-fold 
transformation, there are additional possibilities for observing this property. 

Solutions of the initial problem for f = fi are known [F . Having constructed the transformation 
operator by using ^I7\ and (jSH), we find the solutions Auit) and A22(t) of the set of equations 
after the second step of transformation. Note that the normalization of solutions changes as a result 
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Figure 2: The electric component of the external electromagnetic field after a two-fold transformation 
with the parameter values /o = — 1, C = 0, W21 = 2, Eq = 0.7. 

of the transformation + ^22 — ('C^ + -^i)^- Imposing on the solutions the initial conditions 
^12(0) = 1 and ^22(0) = 0, we determine the probability to populate the excited level after a two- 
fold transformation. The expression for the probability is rather involved and we will not present it 
here. We merely note that this quantity oscillates with time due to the presence of trigonometric 
functions. 

Let us consider the case C = in more detail. It turns out that under the condition 

5/0' - Wf^e + e = 0, (50) 

the oscillations in the time dependence of the probabihty disappear and it acquires a monotonous 
character: 



P2(t) = 1^221' = 

2000t^e^(45(25 + llTs) + 60(123 + 55v^)t^^^ + 32(1525 + 682\/5)tY) 
(5 + v^)5(225 + 540(5 + 2v^)t2^2 + 240(9 + 4:V5)t^^^ + 64(85 + 38^5)^6^6) ' 

Therefore, according to (j^UI) . as distinct from one- fold transformations, for two-fold transformations 

we can indicate two possibihties, /g = ^1 + and /q = ^1 — , which transform the 

probability to populate the excited level from an oscillating function of time to a monotonous one. 
The plots of probability in these cases are given by Figs. 3 and 4 (thick lines). 

We observe another interesting effect by plotting simultaneously the probability to populate the 
excited level before the transformation and after the first and second transformation steps. Fig. 5 
presents the probability plots for a particular value of the parameter /q = which corresponds to 
the case when the probability of population resulting from the one-fold transformation is described 
by a monotonous function jS] 

P.WH^.P^^. (51) 
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10 30 t 

Figure 3: The probability of population of the excited level before the transformation, and after the 
first and second transformations for /q = fl + ^y^j, ^ = 1- 




10 30 



Figure 4: The probability of population of the excited level before the transformation, and after the 
first and second transformations for /q = ^1 — ^5^^, ^ = ^■ 

The thick line presents the plot of probability ()5H) . Its maximal value coincides with the probability 
before the transformation, i.e., with the probability for the usual Rabi oscillations (the dotted line). 
The thin line presents the probability resulting from a chain of two transformations. We see from 
Fig. 5 that the oscillations in this case are in antiphase with the Rabi oscillations. 

Returning to Figs. 3 and 4, note that they also present the plots for /q = ^1 + j and 

/o = (^1 ~ Ts)' respectively. For such values of /o, the probability to populate the excited level 
after the two-fold transformation becomes a monotonous function and for t — 00 the probability 
values after the second transformation (the thick line) coincide with either the maxima or the minima 
of probability oscillations resulting from the one-fold transformation (the thin line). The probability 
before the transformation (the dotted line) oscillates in this case in phase with the probability after 




Figure 5: The probability of population of the excited level before the transformation, and after the 
first and second transformations for /q = ^ = 1. 
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Figure 6: The probability of population of the excited level as a result of a two-fold transformation 
for different factorization constants: /o = 1, -Ri = 0.993, R2 = 0.995, ^ = 0.65. 



the one-fold transformation. 

Omitting the details, let us consider the behavior of the probability to populate the excited 
level obtained as a result of a two-fold transformation with different factorization constants. The 
transformed potential /2 is obtained in 



/2 



fo - 1 [Rl - Rl) 
_ Uk_ _ 1 - iqk 

Uk 1 + iqk ' 



R2 [Bl - {Blf] - R, [Bl - [Blf] 



\R2B2 — RiBi\ 



A; = 1,2. 



(52) 



Using ()42|) and we construct the intertwining operator L02 and then apply it to find solutions 
A12 (t) and A22 (t) of equation ((7j) with potential (jH2I)- In terms of the functions B^, the operator 
Lq2 is written as 



L 



where 



^1 



^2 













S2 








^2 



S2 



Rl 



R2 



R2B2 

^/o + /o + 



RiBl ' 

R1R2 {R2BI 



RiBl) - ifo {R2 - Rl] 



1X2^2 — ^l-^i 



Let 



Rk cos Wkt + hzuic sin w^t — /q 



^k = 2a/ /o 



Rl k 



1,2. 



fo cos Wkt - Rk 

Acting by L02 on solutions of the initial problem ((7j), we obtain its solutions for /o = f2- Imposing 
on these solutions the initial conditions A12 (0) = 1 and A22 (0) = 0, we determine the probability 
to populate the excited level after a two-fold transformation. In general, the behavior of probability 
has an oscillating character. It is important to note that we now observe two kinds of oscillations: 
small-amplitude and large-frequency (fast) oscillations occur at the background of a large amplitude 
and small frequency (slow) oscillations (see Fig. 6). Furthermore, at certain values of parameters, 
the probability for fast oscillations changes within a narrow range (on Fig. 6 between 0.8 and 0.95) 
with the maximal value close to unity. This behavior of probability continues for a long time (on Fig. 
6 for approximately 170 time units that we use). However, the time during which the probability 
of a transition to an exited state is less than 0.5 (in the present case it approximately is equal to 
0.35) is less than the indicated interval by several orders of units (approximately two units in the 
current case). This means that in the course of time the most part of two-level atoms are found in 
the excited state, which may lead to the appearance of an inverse population, as well as to a possible 
appearance of the lasing property of an ensemble of two-level atoms placed in such kind of field. 
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6 Conclusion 



We have considered chains of SUSY transformations for the spin equation, with an apphcation to 
a two-level atom in an external electromagnetic field with a possible time-dependence of the field 
frequency. It is shown that when a chain of n transformations is replaced by a single n-fold transfor- 
mation one can construct a super-Hamiltonian and supercharges which close a polynomial superalge- 
bra. This fact permits us to state that a polynomial pseudo-supersymmetry can be associated with a 
physical system, being a two-level atom in the present case. It is discovered that, as a result of a two- 
fold transformation, a certain type of time behavior of the field frequency implies the disappearance 
of time oscillations of the population probability for the excited level, and the probability becomes a 
monotonously growing function of time with the limiting value which can exceed 1/2. This feature 
permits us to suppose that an ensemble of two-level atoms placed in this specific electromagnetic 
field may acquire an inverse population and exhibit lasing properties. 
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